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1. Introduction 

1.1. WHY DO WE STUDY QUANTUM TRANSPORT? 

As a commonplace, we could say that all simple problems in physics have 
been solved a long time ago, and that we are tempted to turn to the chal- 
lenging field called, rather pompously, "wave transport in complex sys- 
tems". Here, "wave" is ment to emphasize the influence of interference 
(which is truly quantum mechanical for massive particles like electrons or 
atoms). "Transport" implies that we are interested in situations out of ther- 
modynamic equilibrium (but not too far, so that linear response theory ap- 
plies). Finally, a system will be called "complex" whenever it is disordered, 
or strongly interacting, or chaotic. Our choice of this field is motivated by 
two aspects: on the side of applied physics, all remote sensing techniques 
need to incorporate the multiple scattering of waves in turbid media, and 
the miniaturisation process in semi-conductor industry arrives at length 
scales where the control of quantum interference becomes crucial. On the 
academic side of more fundamental physics, we would like to understand 
and enjoy the predictive power of the best physical theory available today: 
quantum theory. 

The classical paradigm of transport in a disordered environment is dif- 
fusion. In general, a quantity n(x, t) in space and time (think of a particle 
concentration) at equilibrium is distributed with a constant density uq. If 
the quantity is locally conserved, a small variation 6n{x, t) then obeys the 
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diffusion equation, [dt — Do'V'^)6n{x, t) = 0; Dq is tlie diffusion constant. 
Its solution in Fourier space 

6h{q, t) = 6h{q, 0) exp(-L>og^t) (1) 

shows tliat very smootli fluctuations in real space (with small wave numbers 
q ^ 0) are removed on very long time scales l/Dgq^ — > oo. Indeed, because 
of the local conservation law, a small surplus 6n > cannot just simply 
disappear, but has to reach a place far away with a depleted density Sn < 
in order to restore the equilibrium situation Sn = 0. 

Classical Boltzmann diffusion theory has been used successfully to de- 
scribe electric conduction in metals (first by Drude), or the diffusion of 
light intensity through stellar atmospheres or interstellar clouds (notably 
by Schuster and Chandrasckhar) [1]. However, at the microscopic level, one 
has to deal with waves. As an example, the quantum picture of electrons 
moving in a perfectly periodic crystal pinpoints the key role of interference: 
an electron, initially confined in a well, can resonantly tunnel through the 
lattice, yielding Bloch energy bands. In trying to understand the interplay 
between disorder and interference, P. Anderson showed that sufficiently 
strong disorder can suppress the quantum diffusion (leading to a vanishing 
diffusion constant D = 0), a phenomenom baptised localisation [2]. Later, 
it was realized that even far from the regime of (strong) localisation, dif- 
fusive transport is affected by interference: this so-called weak localisation 
reduces the Boltzmann diffusion constant, D = Dq—SD, by the constructive 
pairwise interference of amplitudes associated with time-reversed scattering 
paths. Mesoscopic physics, namely the study of interference effects in wave 
transport through random media, was born [3]. 

1.2. WHY DO WE USE PHOTONS AND ATOMS? 

The theory of localisation has been first developed in the condensed mat- 
ter physics community. But as electrons are charged particles with a very 
strong and long-range Coulomb interaction, the pure one-particle picture of 
Anderson localisation has never been observed experimentally as far as we 
know. On the other hand, the radar physics community first had realized 
that interference of counter-propagating amplitudes can play an impor- 
tant role in the multiple scattering of electromagnetic waves [4] . Taking the 
better of the two worlds, S. John and P. Anderson suggested to study local- 
isation using light or other non-interacting classical waves. Light scattering 
allows the use of modern lasers with excellent coherence properties as well 
as an accurate analysis of direction and polarisation. A large number of 
turbid media has thus been studied, from Saturn's rings to semi-conductor 
powders [5]. 
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Figure 1. Experimental CBS enhancement a vs. the angle 6 with respect to the backscat- 
tering direction for a classical medium (Styrofoam, left) and cold rubidium atoms (right), 
in the channel of preserved helicity {h\\h, top) and flipped helicity (h ± h, bottom). Data 
generously provided by G. Labeyrie. 



Today, laser-cooled atoms as light scatterers come close to a theoreti- 
cian's dream team: they are perfectly identical (monodisperse) point parti- 
cles. But more than that (theoretical convenience hardly ever justifies ex- 
pensive experiments): their specific properties permit to study new regimes 
which are characterised, for example, by their quantum internal structure 
(as discussed in the present contribution), the finite width of the atomic 
fluorescence spectrum, non-linearities such as the saturation of an atomic 
transition, and the mechanical acceleration of the atom due to light scat- 
tering. Even more, at sufficiently low temperature, the atoms themselves 
become matter waves and can in turn be used to probe quantum transport, 
all the way down to the degenerate quantum regime of Bose-Einstein con- 
densation (see, for example, the recent beautiful experimental realisation 
of the Mott-Hubbard quantum phase transition by the Munich BEC group 
[6]). 

1.3. WHY DO WE NEED TO DESCRIBE AN INTERNAL STRUCTURE? 

The first experiment of coherent backscattering of light from a cloud of cold 
atoms in 1999 [7, 8] yielded a surprising result (see Fig. 1): the observed 
interference peak in the channel of preserved helicity [h\\h) shows only a 
maximum enhancement of about 1.05, far below the classically expected 
factor of 2.0 due to reciprocity [9]. Also, this h\\h enhancement is much 
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Figure 2. Left: a cloud of cold atoms as point scatterers with fixed random positions 
r and randomly oriented total angular momenta J. Right: a zoom into the energy level 
scheme of a resonant degenerate dipole transition, here the example of J = 1, Je = 2. 
(5 = — cjo is the detuning of the probe light from the atomic resonance frequency tJo , and 
r is the natural width of the excited atomic state. The polarisation of scattered polarised 
photons (full and dotted lines) is coupled to the internal magnetic quantum numbers m. 

smaller than the measured value of 1.2 in the channel of flipped helicity 
(/i_L/i). It was soon realized that the degeneracy of the probed atomic 
dipole transition J = 3 — > Je = 4 is responsible for an imbalance of CBS 
amplitudes and therefore reduces the measurable enhancement factor [10]. 
We thus have to generalise the theory of the multiple scattering of polarised 
light by point dipoles to the case of an arbitrary atomic transition J 
Je [11]. This theory indeed explains the observed enhancement factors and 
shall be described in the following. 

2. Multiple scattering of a photon by atoms with internal degen- 
eracy 

2.1. THE ONE-PHOTON TRANSITION MATRIX 

Consider a cloud of laser-cooled atoms confined in a standard magneto- 
optical trap. The cooling is such that their velocity spread v is much smaller 
than the Doppler velocity F/fc (T is the natural width of the excited atomic 
state). Therefore, we can neglect the Doppler effect and may assume that 
the atoms' positions r^, a = 1, . . . ,N, remain fixed on the light-scattering 
time scale. On the other hand, the velocity spread should be much larger 
than the recoil velocity Vj-^c = hk/M (where M is the atomic mass) for the 
scattering of a photon of wave-vector k. This allows us to treat the positions 
as classical random variables and to follow the standard diagrammatic ap- 
proach to describe multiple scattering (see [12] and references therein). The 
CBS probe beam with incident wave- vector k, polarisation s and frequency 
iv excites a closed atomic dipole transition defined by a ground state with 
total angular momentum J and an excited state Jg with frequency ujq. In 
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the absence of a magnetic field, these two levels with internal quantum 
numbers m and me are respectively (2J+ 1)- and (2Je + l)-fold degenerate 
(see Fig. 2). 

In order to minimise saturation effects and optical pumping, the exper- 
iment is performed at low laser intensity, so that we can consider the light 
scattering in the limit of one-photon Fock states \k£) (in this notation, the 
transversality e • = is understood). The transition amplitude for the 
scattering of an incident photon \k£) into an emitted photon {k'e') by a 
single atom is 

k £ -> Jm' 

^n,,^^p^^*C^ ={Jm'\e' ■t{io)-e\Jm)e'^''-'''>'-. (2) 

Note that the exponential with the classical external degrees of freedom r 
is factorised from the matrix element with the internal quantum numbers. 
Therefore, the usual multiple scattering formalism applies [12]. But we have 
to analyse carefully the role of the internal degrees of freedom. We see 
in (2) that the incident and emitted polarisation vectors e and e' (the 
bar denotes complex conjugation) are coupled by the transition matrix 
tij{m,m';ijj) = {Jm\tij{uj)\Jm). As any 3x3 matrix, it can be decomposed 
into its irreducible components with respect to rotations, 



its scalar part or trace t^^\ its antisymmetric part t^^^ and its traceless 
symmetric part t^'^^ (summation over repeated indices is understood). It 
is easy to verify that in the non-degenerate limit J = 0, corresponding 
to the classical case of an isotropic dipolc, only the scalar component tf^^ 
exists. So our task in the following will be to determine how the non-scalar 

(1 2) 

components t-' influence the light propagation. 

2.2. AVERAGE PROPAGATION INSIDE THE EFFECTIVE MEDIUM 

In a disordered environment, the actual propagation depends on the pre- 
cise realisation of disorder and yields, for example, an interference "finger- 
print" or speckle-pattern. Universal properties can only be obtained on 
average. The average propagation of a light mode \k£) inside the medium 
is described by the average propagator {Gij{k,oj)) where the average (. . . ) 
traces out the matter degrees of freedom. Following the standard procedure, 
we determine the average propagator in terms of the vacuum propagator 
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go and the self-energy S with the aid of the Dyson equation (tensor indices 
are omitted for brevity) 

(G)=5o + 5oS(G) = ^^. (4) 

For a sufficiently dilute scattering medium such that nX^ <C 1 (n is the 
atomic number density with typical values of 6 x 10^'^ cm~^, and A is the op- 
tical wavelength, typically 780 nm in the experiments done so far), the inde- 
pendent scattering approximation applies, and the self-energy is simply the 
average single-scattering t-matrix: T,ij{uj) = N {tij{uj)). We suppose that 
the atomic scatterers are distributed uniformly over their internal states 
such that rotational invariance is restored on average: Sij(w) = 6ij'E{ijj). 
This is intuitively clear since under a scalar average only the scalar com- 
ponent can survive. Now, it can be checked that the scalar component is 
precisely the one that mimics a classical isotropic dipole: 

Here, the quantum internal structure only enters through the scalar factor 

Mj = (2Jc + 1)/3(2J + 1) with the non-degenerate limit Mq = 1. The 
elastic mean free time is defined as r = — (2ImS(a;))~^ and yields the mean 
free path £ = r in our units, where c = 1. By virtue of the optical theorem, 
the mean free path is related to the total cross-section for elastic scattering 
by the usual Boltzmann expression i = 1/na. 

Obviously, the quantum internal structure (or contribution of non-scalar 
components of the t-matrix) has disappeared under the scalar average over 
internal states, and we recover a standard scalar theory. Have we been too 
optimistic in hoping that the internal structure can explain the dramatic 
reduction of interference in the CBS signal? Well, of course not: the experi- 
mental signal in Fig. 1 is the average intensity which must be distinguished 
from the square of the average amplitude. More technically speaking, af- 
ter calculating the average amplitude {G) with the Dyson equation (4), we 
have to turn to the average intensity (GG). 

2.3. THE AVERAGE INTENSITY VERTEX 

In a dilute medium <C 1, the building block for the multiple scatter- 
ing series is the average single scattering intensity vertex {tijiki) [13]. It 
connects two amplitudes with their respective polarisation vectors and can 
thus be written as a rank-four tensor, 

{tij{u)iki{u)) = u{oo) I I . (6) 
l-^k 
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Here, the four point polarisation vertex can be calculated analytically using 
standard methods of irreducible tensor operators [14] . It is simply given as 
the sum of the three pairwise contractions, 

^ 11 ^^ = wiSijdki + W2Sikdji + wsdiiSjk (7) 

where the coefficients wi = -^{sq — S2), W2 = ^{s2 — si) and W3 = ^{s2 + si) 
are given in terms of the squared 6 J-symbols 

.K = 3(2Je + l)| J J j^f (8) 

associated with the irreducible t-matrix components t^^^ of order K = 
0, 1, 2. For the classical isotropic dipole J = 0, the coefficients become 
{wi,W2, ws) = (1, 0, 0) and yield the purely horizontal contraction 5ij6ki- We 
see therefore that including the quantum internal structure is equivalent to 
replacing the simple dotted line in the usual diagrams [12] by the generalised 
vertex (7) with the richer topology of a ribbon: 




2.4. SUMMATION OF LADDER AND CROSSED SERIES 



Following the standard diagrammatic approach, we have to sum the so- 
called series of ladder diagrams 



L 



+ ! ! + ! ! + 

(g^^g) (g^^gM^g) 



(10) 



that describe the average intensity neglecting all interference terms: the 
direct amplitude (upper line) is scattered by the same scatterers as the 
conjugate amplitude (lower line), as indicated by the dotted lines. At least 
formally, one recognises a geometrical series that may be summed up ana- 
lytically once the single scattering vertex (first term on the r.h.s.) and the 
square of the average propagators (thick lines) are known. The interference 
correction associated with CBS and weak localisation is contained in the 
so-called maximally crossed diagrams 



C 



(g^^ g^^ 

'■^<" + + 
(g— — ig) (giii__(g^_i<g 



(11) 
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that describe the propagation of the direct and conjugate amphtude along 
the same scattering path, but in opposite directions. For classical point scat- 
terers, the crossed and ladder diagrams are closely related by reciprocity: 
by turning around the lower line of a maximally crossed diagram, the con- 
necting lines straighten out and yield the corresponding ladder diagram. By 
carefully counting incident and emitted momenta and polarisation indices, 
one shows that the diagrams are rigorously equal for scattering in the back- 
ward direction (k' = — fe) and in the parallel polarisation channels e' = e. 
This identity justifies the CBS enhancement factor of 2.0 in the helicity- 
preserving polarisation channel where, for spherically-symmetric scatterers, 
the single scattering contribution is absent. 

In our case of a quantum internal structure, see (9) , the classical vertex 
has to be substituted by the ribbon vertex. But now the correspondence 
between ladder and crossed diagrams is spoilt: when the bottom line of a 
crossed diagram is turned around, the connecting ribbons are twisted: 



The evident difference between the twisted and the straight ribbon shows 
that the quantum internal structure of the atomic scatterer indeed affects 
the interference corrections to the average intensity. More quantitatively, 
twisting the vertex is equivalent to the exchange of the coefficients W2 and 
u)3 associated with the diagonal and vertical contractions. Once we have 
calculated the ladder series, we can then obtain the crossed series by simply 
replacing W2 ^ (and rearranging the momenta and polarisation vectors 
as in the classical case). 

In order to sum the geometrical ladder series, we have to determine the 
eigenvalues of the atomic vertex with respect to the "horizontal" direction 
of summation. Using a basis of projectors T*-^) onto irreducible eigenmodes 
[eje/]^^-* of the field polarisation matrix, we obtain a decomposition of the 




(12) 



form 




K 



(13) 



where the eigenvalues \k are simple functions of the Wi'. 



Ao = l, 



Xl=Wi- W2, 



\2 = U!l + W2. 



(14) 
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Taking advantage of our simple substitution rule W2 ^ w^, the twisted 
vertex of the crossed series then gives crossed eigenvalues XK such that 

= T.XK{J,Je)TZl (15) 
K 

Having diagonalised the scattering vertex, we have to treat also the trans- 
verse propagation between atoms. The actual calculation is rather cumber- 
some because it involves momentum-dependent eigenvalues and projectors; 
details can be found in [15]. But our approach permits us to sum analyti- 
cally the ladder and crossed series, for the full transverse vector field with 
arbitrary polarisation, and for arbitrary atomic transitions. It is a first step 
towards a generalisation of the existing multiple scattering theories, either 
of scalar waves by anisotropic point-scatterers [16, 17], or of vector waves 
by isotropic Rayleigh scatterers [18, 19]. 




3. Bulk transport properties 

3.1. DIFFUSION AND DEPOLARISATION 

The summed ladder propagator describes the average intensity distribu- 
tion inside the bulk medium (i.e. in the absence of boundaries). To gain 
qualitative insight, we can simplify the exact expressions by the diffusion 
approximation (retaining terms of order q^). The crucial ingredients are the 
atomic eigenvalues Xk and xk as well as the eigenvalues bo = 1, bi = ^ and 
b2 = jQ associated with the transverse propagation. The ladder propagator 
in the long-time limit then takes the form 

L{q, t)=Y,lK exp [-DkqH - t/T^{K)\ . (16) 

K 

Apart from the factor Ik with no importance here, we first notice an ex- 
ponential decay that was anticipated above in the general diffusion picture 
(1). We obtain the diffusion constant Dk ~ Dq = in terms of the 

well-known transport velocity utr = ^/Ttr of resonant point scatterers [12]. 
But there is a second exponential decay, with a characteristic polarisation 
relaxation time 

The scalar field mode {K = 0) describes the total intensity. Its relaxation 
time is infinite (Tp(0) = oo) for all atomic transitions (since 6o = -^o = 1)) 
and we recover a purely diffusive behaviour as required by energy conser- 
vation. The non-scalar field modes are exponentially damped on finite, and 
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Figure 3. The transversality of propagation depolarises the multiply scattered light. This 
depolarisation is enhanced by an atomic quantum internal structure J > (as discussed 
in sec. 3.1). Furthermore, the internal degeneracy leads to dephasing of the interference 
correction (as discussed in sec. 3.2). 



in fact rather short time scales rp(l,2) < rtr: this simply says that a well- 
defined field polarisation is scrambled in the course of multiple scattering. 
With increasing degeneracy J > of the atomic transition, these times 
get even shorter (A;^ < 1), which nicely confirms the intuitive picture that 
random transitions between different atomic Zeeman substates enhance the 
depolarisation. 

3.2. WEAK LOCALISATION AND STRONG DEPHASING 

Under the same approximations, the sum of crossed diagrams yields a con- 
tribution 

C{qc, t) = ^CK exp [-DKq'^ct - t/T^{K) - t/T^{K)\ (18) 
K 

as a function of the total momentum qc = \\k + k'\\. One first recognises 
the same exponential decay as in the ladder contribution (17), indicating 
that depolarisation also affects the coherent contribution. But there is an 
additional source of exponential damping described by dephasing times [20] 

which we define precisely as the damping times with respect to the inco- 
herent depolarisation times. Of particular interest is the dephasing time 
T0(O) of the scalar mode or intensity. For atomic transitions of the type 
Je = J + 1, we find explicitly 

The interference is only preserved (r^ = oo) for classical dipoles J = 0. 
For the least possible degeneracy J = i, the dephasing time is already as 
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short as rtr/2 and decreases as 1/J^. Even if the exact expression becomes 
meaningless for larger J (since the diffusion approximation is certainly not 
justified on time scales shorter than the transport time), it is evident that 
the internal structure destroys the interference very effectively. Therefore, 
in the presence of a quantum internal structure, we only expect a (very, 
very) weak localisation correction to the Boltzmann diffusion constant of 
light propagation. 

4. Coherent backscattering 

We wish to calculate the CBS peak [21] analytically for arbitrary atomic 
transitions and therefore choose the simplest possible geometry of the scat- 
tering medium, a semi-infinite half-space. Having calculated the bulk prop- 
agator FqIvi — r2) (either the ladder or the crossed contribution), we define 
the corresponding propagator for the semi-infinite half space by using the 
method of images, F{ri,r2) = -Fo(^i — ^2) — Fo{fi — f2')- The image point 
2' is defined with respect to a mirror plane lying at a distance zq outside 
the boundary of the medium. Following the habits of the field, we use the 
so-called "skin layer depth" zq 0.7121 £ pertaining to the exact solu- 
tion of the homogeneous Milne equation for vector waves and point dipole 
scatterers [18, 19]. 

By this approach, we can calculate the CBS enhancement factors and 
peak shapes beyond the diffusion approximation (which is crucial whenever 
only short paths contribute). The CBS enhancement factor a is plotted in 
Fig. 4 as a function of J for transitions of the type J J + 1. In the 
classical limit J = 0, we recover values of a = 2.00, 1.76, 1.24 and 1.12 for 
the different polarisation channels which are in excellent agreement with 
the exact values obtained by a solution of the vector Milne equation by the 
Wiener-Hopf method [18, 19]. This indicates that the method of images can 
be used with success even for signals involving short scattering paths (like 
in the perpendicular channels) provided that the exact propagator (beyond 
the diffusion approximation) be used. 

Figure 4 shows that the least internal degeneracy reduces the CBS inter- 
ference dramatically: the perfect factor of 2.0 in the h\\h channel plunges 
down to 1.04, well below the other three polarisation channels. At this 
point, our theory indeed explains the astonishing experimental result shown 
in Fig. 1 that has motivated this work. Furthermore, an experiment using 
cold Strontium atoms without internal degeneracy (J = 0) has recently 
confirmed that excellent enhancement factors in agreement with the theo- 
retical predictions are obtained [22]. 

Figure 5 shows a comparison between calculated CBS peak shapes for 
isotropic point scatterers (left side) and atoms with J = 3 ^ Je = 4 
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Figure 4- Calculated CBS enhancement factors at exact backscattering from a homo- 
geneous half-space of atoms with a degenerate dipole transition J —f J + 1 for the four 
usual channels of preserved (||) and flipped (_L) linear polarisation (l) or helicity (h). 
Contrary to the classical case, for atoms with J > the best enhancement is expected 
in the hJ-h channel. 



(right side). The full lines are the sum of all scattering orders, and the 
dotted lines indicate the double scattering contribution (which is known in 
closed form [14]). Whereas the CBS peak in the classical case and parallel 
polarisation channels contains very high orders of scattering (corresponding 
to long scattering paths), the atomic internal degeneracy cuts off these 
contributions (as indicated already by the dephasing times (19)) and yields 
only very small peaks. 

A quantitative comparison of the theory to the experimental results 
needs to take into account the finite geometry of the actual atomic cloud 
(roughly spherical, with a Gaussian density distribution). This means that 
analytical results are out of reach, and a numerical approach has to be 
taken. A Monte Carlo simulation of photon trajectories in various geome- 
tries has been realized by D. Delande and yields results which are in good 
agreement with the experimental data [23]. 
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J = 0, Je = 1: J = 3, = 4: 




-10 -5 5 10 -10 -5 5 10 

MO kie 

Figure 5. Calculated CBS enhancement factor a as a function of the reduced backscat- 
tering angle ki9 for isotropic dipole scatterer (left) and atoms with a degenerate 
J = 3 — > Je = 4 transition (right). Full line: sum of all scattering orders. Dotted line: 
double scattering contribution. 

5. Concluding remarks and acknowledgements 

In summary, we have presented an analytical theory of the multiple scat- 
tering of polarised photons by resonant atomic dipole transitions with ar- 
bitrary degeneracy. We have shown how the usual diagrammatic approach 
can be generalised by using an intensity vertex with a "ribbon" topology 
that breaks the equivalence of ladder and crossed diagrams. The theoretical 
CBS peak heights reproduce the experimental results: the quantum inter- 
nal structure indeed reduces the CBS interference drastically. Inasmuch as 
weak localisation acts as a precursor for Anderson localisation, our results 
indicate that in order to reach the strong localisation regime with cold 
atoms, the use of a non-degenerate transition is highly recommendable. 
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We hope to have demonstrated that atoms are Ught scatterers with in- 
triguing interference properties. They permit to set up a microscopic theory 
for diffusion and weak locahsation and thus promise to be a source of fur- 
ther inspiration for both fields of "disordered systems/multiple scattering" 
and "atomic physics/quantum optics". 

It is a pleasure to acknowledge the invaluable help by Dominique De- 
landc (who is to be credited for the "exact-image" method), many inspiring 
discussions with Eric Akkermans (who is to be credited for the "dephas- 
ing" interpretation), and a critical reading of the manuscript by Andreas 
Buchleitner. 
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